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1.  Introduction 


This  paper  examines  some  aspects  of  a stochastic  process  which  might 
be  vised  to  model  the  mobility  of  users  in  a mobile  communication  network. 
Specifically,  we  model  the  notion  of  a user  as  a two-dimensional  sym- 
metric random  walk.  This  might  be  eui  appropriate  model  for  a patrol  car 
in  a city.  We  assume  that  there  is  a controller  who  needs  to  be  kept  in- 
formed of  the  locations  of  the  various  users  to  within  some  tolerance.  Here 
we  will  require  that  each  user  notify  the  controller  of  its  exact  location 
whenever  it  hits  a square  boundary  of  a certain  size  centered  on  the  point 
it  occupied  (Point  C in  figure  1)  the  last  time  it  communicated  its  location 
to  the  controller. 


User  informs  controller 
of  its  location  os  soon 
os  it  hits  this  square 
boundary 


Figure  1 


Since  the  controller  knows  point  C,  the  user  need  only  transmit  an 
index  which  indicates  the  point  at  which  it  hits  the  square  boundary.  In 
the  following  sections,  we  will  derive  expression  for  (1)  the  expected 
number  of  steps  or  time  a user  takes  to  reach  the  boundary  starting  from 
the  center  C and  (2)  the  probabilities  with  which  the  user  hits  the  dif- 
ferent point  on  the  boundary  and  the  corresponding  entropy.  It  is  clear 
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that  the  smaller  the  size  of  the  square  boundary,  the  better  informed 
is  the  controller  of  the  location  of  the  user.  On  the  other  hand,  one 
would  surmise  that  the  user  will  have  to  transmit  more  data  back.  This 
trade-off  is  discussed  in  section  IV. 

II . Expected  time  between  user  transmissions. 

In  this  section,  we  derive  an  expression  for  the  expected  time  for  a 

two-dimensional  symmetric  random  walk  starting  from  the  center  of  a square 

boundary  to  reach  some  point  on  the  boundary.  For  convenience,  we  consider 

2 

the  square  to  have  length  it  emd  subdivide  it  into  n small  squares  as  shown 
in  figure  2 for  the  case  n=4.  We  assume  n to  be  an  even  integer. 
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Figure  2 


Let  us  denote  the  e3q)ected  time  to  hit  the  boundary  starting  from  the  point 
(x,y)  by  t 


. Of  course  x and  y take  on  values  in  the  set  {O,  — , — , ...f- 
x,y  ^ n ' n n 


We  are  primarily  interested  in  though  it  will  be  f2u.rly  easy  to 

write  ejq>ressions  for  t for  arbitrary  x and  y. 

x,y 

starting  from  (x,y) , the  user  can  move  to  (x+h,y) , (x-h,y) , (x,y+h)  or 
{x,y-h)  where  h = 7T/n  and  each  motion  occurs  with  equal  probability  1/4.  The 
conditional  expected  time  to  hit  the  boundary  assuming  the  first  step  takes  the 
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user  to  (x',y')  is  t^,  This  argument  shows  that  t^  ^ satisfies  the 

difference  equation 


t =Tt,  +Tt  u+Tt  ^^-  + 1- 

x,y  4 x-h,y  4 x+h,y  4 x,y-h  4 x,y+h 


where  0<x<Ti,  0<y<TT  with  the  boundary  conditions 


(la) 


t„  =t  =t„=t  =0. 

0»y  Tr,y  x,0  x,Ti 


In  order  to  solve  this  system  of  equations,  we  first  consider  the 
eigenvalue  problem 


(lb) 


4Xt  =t.  +t^.  +t  ,+t 

x,y  x-h,y  x+h,y  x,y-h  x,y+h 


(3) 


associated  with  the  homogeneous  part  of  (1) . It  should  be  observed  that  as 

, ..  . 2tt  (n-l)~i,_.  . , ,,2..,  ,,2 

X and  y range  over  the  set  i— , — , H3)  represents  an  (n-1)  x(n-l) 

n n n 

matrix  eigenvalue  problem.  The  eigenvalues  and  eigenvectors  of  this  system 
are  given  in  (1,  p.289].  The  eigenvalues  are 


X _ = •:7  (cos  ph  + cos  qh)  , p,q  = 1,2,..., n-1 

p»q  2 


and  the  corresponding  eigenvectors  are 


(4) 


p»q 


1,1 

^1,2 


“l,n-l 

'^2  ,1 


2, n-1 


L n-l,n-lj 
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where  u = sin  prh  sin  qsh,  p,q  = This  can  be  verified  by 

i ^ s 


setting  t =sin  px  sin  qy  and  noting  that 


2 sin  px  cos  ph  = sin  p(x+h)  + sin  p(x-h) . 


(6) 


We  now  show  that  the  eigenvectors  defined  by  (5)  are  mutually  orthogonal, 
that  is. 


n-1 


D = I sin  prh  sin  qsh  sin  p'rh  sin  q'sh  = 0 
r,s=l 


(7) 


if  p'  j/i  p or  q'  5<  q. 


n-1  ^ n-1 

D = y sin  pr  — sin  p'r  — . J 
r=l  s=l 


ir  . , IT 

sin  qs  — sin  q's  — 


(8) 


Thus  to  show  that  the  eigenvectors  are  orthogonal  it  suffices  to  show  that 
if  P'  7^  Pf 


n-1 


y sin  pr  — sin  p'r  — = 0. 
, n n 

r=l 


(9) 


n-1  n-1 

y sin  pr  — sin  p'r  — = 7 zr  [cos(p-p')r  — - cos(p+p')r  — ] (10) 

**  1^  n ^ X T1  M 


••  ‘'i  2 

r=l  r=l 


From  [2,  p.78,  Eq.  (418)], 


n-1 


7 cos  r(p-p')  — = cos  (p-p’)  ? sin  (n-1) (p-p')  cosec  (p-p')  ~ 
, n ^ zr 


r»l 


2n 

(11) 


'COS  (p-p')  — [ sin  (p-p')  — cos  (p-p')  - cos  (p-p*)  — sin  (p-p*)  — ] 


2n' 


cosec  (p-p*)  j. 


(12) 


2 / IX 

» - cos  (p-p  ) ~ . 


n-1 

r TT 

I COS  r (p+p' ) — = 

3:=!  ”■ 


2,  . TT 
- cos  (p+p')  J 


Finally,  substituting  (13)  and  (14)  in  (10)  yields 


I sin  pr  ^ sin  p'r  ^ j [ -cos^(p-p')  j + cos^(p+p')  j 1 (15) 


= 7 [ {-1-  cos(p-p’)Tr}  + {l  + cos(p+p’)7r}  ] 


— I - 2 sin  pir  sin  p'TT]  = 0 


It  is  now  fairly  easy  to  obtain  the  solution  of  (1)  in  terms  of  the 

2 

eigenvectors  U • Let  T denote  the  (n-1)  dimensional  vector 

p»q 


^,h 

Hi,2h 

(n-l)h 


2h, (n-l)h 


'(n-l)h,  (n-l)h 


Since  the  ^ are  linearly  independent,  we  can  write 


n-1 

I 6 u 

p.q=l 


T » 


(18) 


I 


L_ 
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for  some  constants  3 to  be  determined.  Similarly  we  can  write  the  (n-1) ' 


vector  1,  all  of  whose  components  are  ones,  as 


n-1 

I 


p,q=l 


a u 
p»q  p»q 


(19) 


Taking  the  dot  product  of  both  sides  of  (19)  with  U.  . and  recalling  the 


fact  that  the  eigenvectors  U are  orthogonal,  we  obtain 

p»q 


1.  u ^ = a . u .| 


(20) 


2 2 2 
where  |u.  -1  = I sin  irh  sin  jsh. 

r,s=l 


(21) 


2 2 

It  Ccin  be  verified  [2,  p.82]  that  |u.  .(  = n /4.  Also, 

» D 


n-1 


1.  U.  . = Y 3®^ 

r,s=l 


(22) 


= -i-  [l-(-l)^]  [l-(-l)^]  cot  g cot  g 


(23) 


In  (23)  we  have  used  the  fact  [2,  p.78]  that 


n-1 


r=l 


sin  irh  = [l-(-l)^]  cot  ^ . 

2 2n 


(24) 


Thus  from  (20)  we  can  write 


[l-(-l)^]  [l-(-l)^]  cot  ^ cot  ^ 


i.j 


il 


(25) 
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Now,  if  we  denote  the  matrix  associated  with  the  homogeneous  part  of 
(1)  by  A,  then  solving  (1)  is  equivalent  to  solving  AT  = IT  - 'J,  i.e. 


BT 


(26) 


2 2 

where  B = A - I and  I is  the  (n-1)  x(n-l)  identity  matrix.  Since 


BU  = AU  - U = (X  -1)U 
p,q  p,q  p,q  p,q  p,q 


(27) 


the  eigenvalues  of  B are  (X  -1)  and  its  eigenvectors  are  U . From  (26) 

p / q p ' ^ 


and  (18)  we  obtain 


n-1 


BT 

Therefore 


R (X  -1)U  = -1  = y - a U 

n.rt  r^_rT  r^rr  i^.rrrv.rT 


p,q  p,q  p,q 


,q=l  P''! 


(28) 


I ' 


M 


I 
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Equation  (31)  gives  us  an  expression  for  the  expected  number  of  steps  a 
symmetric  two-dimensioncil  random  walk  tcikes  to  hit  a square  boundary  of 

* 

size  n starting  from  some  point  inside  the  boundary.  Let  us  denote  by  t^ 
the  expected  time  ■n/2  boundary  starting  from  the  center. 

Then 


* 

t = 
n 


n-1 

I 

p,q=l 


-8  cot  5—  cot  5“  sin  ^ sin  ^ 
2n  2n 2 ^ 

2,  pTT  OTT 

n (cos  + cos  •'* 2) 

n n 


(32) 


Figure  3 shows  t^  for  values  of  n ranging  from  2 to  20. 


Ill . Probabilities  of  boundary  points 

A 

Let  p denote  the  probability  that  the  random  walk  will  first  hit 

* ♦ 

the  square  boundary  at  point  (x  ,Y  ) given  that  it  starts  from  the  point 

* 

(x,y)  . In  this  section  we  will  derive  an  expression  for  tt/2' 

it 

an  argument  similar  to  that  leading  to  (1) , we  see  that  p satisfies  the 
difference  equation 


1*  1*  +i* 

^x,y  4 ^x-h,y  ^ 4 ^x+h,y  ^ 4 ^x,y-h  4 ^x,y+h 


where  0<x<tt,  0<y<7r 


(33a) 


with  the  boundary  conditions 


* * * * * 

p„  =p  =P„=P  =0  except  that  p . . = 1. 

*^0,y  *^Tr,y  ^x,0  ^x,ti  ^ x*,y* 


(33b) 


J 


For  notational  convenience,  let  us  denote  the  point  inside  the  boundary 

2 

immediately  adjacent  to  (x*,  y*)  by  ix' , y').  Also  let  P denote  the  (n-1) 
dimensional  vector 


^h, (n-l)h 


(n-l)h,  (n-Dh 


Then  we  can  write  (33)  as 


AP  = IP  - -r  e . I 

4 x',y' 


where  A is  the  same  matrix  as  in  the  previous  section  and  e , .is  the 

x’  ,y' 

2 

(n-1)  dimensional  vector  with  1 in  the  position  corresponding  to  x' ,y' 

and  O's  elsewhere.  Because  of  the  symmetry,  there  is  no  loss  of  generality 

in  assuming  that  x*  = it  and  ^ y*  £ it  . In  this  case,  x’  = (n-1)  y'  = y* 

eind  the  non-zero  entry  in  e , , is  the  [(n-2)  (n-1)  + X,]-th  entry  where 

X ,y 

Proceeding  in  a manner  completely  analogous  to  the  method  of  the 
previous  section,  and  setting 


we  find  that 


a.  . = 


sin  [i(n-l)h]  sin  jilh 
2 


Also  letting 


n-1 

p = y B U 

p,q=l 


we  find  that  6 


P/q  _ _ 2 sin[p(n-l)h]  sin  qllh 


p.q 


X -1 


p»q 


n (cos  ph  + cos  qh  -2) 


Thus , 


P = I 


■■  . / -.  » ■n'  ■ n T! 

n-1  2 sin  p(n-l)  — sin  qil  — 

” n n 


1 2 _ , TT  . 7r  , , 

p,q=l  n [2  “ (cos  p—  + cos  q ““  ) ] 

n n 


DTT  . qTT 
Sin  • — sin 

n n 

P'TT  . q27T 

Sin  sin  ^ — 


sin  ^ sin  q(n-l) 
n 

. p27T  . qTT 

Sin  * sin 

n n 


p{n-l)7r  . a 

Sin  sin 

n 


From  (39)  we  can  deduce  that 


TT  7T  II 

. n-1  2 sin  p(n-l)  — sin  qZ.  — sin  ^ sin  _ 

” " n n 2 2 


21  31 


•7r/2,Tr/2 


I 

P'‘^=^  n2[2  -(cos  cos  2^)] 


n-l)7T 

n 


t=l  c 
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l 
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Using  (40) , the  entropy  H(n)  in  bits  associated  v;ith  the  random  variable 

indicating  the  boundary  point  the  random  walk  first  hits  was  plotted  as  a 

function  of  n in  figure  4.  For  a square  boundary  of  size  n,  there  are 

4(n-l)  boundary  points  which  could  first  be  hit.  If  we  assume  that  all  these 

points  are  equi-probable,  the  resulting  entropy  H (n)  = 2 + log  (n-1) 

is  clearly  an  upperbound  on  H(n).  It  was  found  numerically  that  for  even 

positive  integers  n less  than  or  equal  to  20,  the  difference  between  H (n) 

max 

and  H(n)  was  monotomically  increasing  but  did  not  exceed  0.18.  The  maximum 
percentage  error  was  below  3%. 

IV.  Discussion 

In  order  to  give  a rough  indication^of  the  trade-off  between  the 
size,  n,  of  the  square  and  the  amovmt  of  information  to  be  transmitted 
back  to  the  controller  every  time  unit,  a plot  of  H(n)/t*  is  shown  in 
figure  5.  From  a practical  viewpoint,  the  trade-off  between  [2  + log2 (n-1) )/t* 
and  n should  be  excimined. 

Finally,  we  note  that  the  symmetric  model  assxamed  here  may  be  refined 
in  many  ways.  For  example,  a more  appropriate  model  of  a patrol  car  might 
exclude  the  possibility  of  the  car  making  a U-tum.  Unfortunately,  such  models 
seem  to  be  hard  to  analyze. 

t Note  that  we  are  neglecting  the  fact  that  there  is  dependence  between 
the  time  of  hitting  the  boundary  and  the  boundary  point  which  is  first 
hit.  A more  accurate  analysis  snould  take  this  into  account. 


